Minigap in a long disordered SNS junction: analytical results 
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We review and refine analytical results on the density of states in a long disordered super- 
conductor-normal-metal-superconductor junction with transparent interfaces. Our analysis includes 
the behavior of the minigap near phase differences zero and tt across the junction, as well as the den- 
sity of states at energies much larger than the minigap but much smaller than the superconducting 
gap. 



A superconductor in contact with a normal metal in- 
duces pairing correlations in the metal, a phenomenon 
known as the proximity effect. One of the most remark- 
able consequences of such induced correlations is the 
appearance of a gap, usually referred to as the "mini- 
gap", in the electronic excitation spectrum of the nor- 
mal metal |^]. A very common setup exhibiting a mini- 
gap is the superconductor-normal-metal-superconductor 
(SNS) junction made of two superconducting leads con- 
nected via a disordered normal layer. The gap in this 
junction is determined by the diffusion time across the 
normal layer and is sensitive to the phase difference 
across the junction, reaching a maximum at zero phase 
difference and vanishing at the phase difference tt 

The appearance of a minigap and its phase dependence 
is well understood in the quasiclassical description. Pro- 
vided the scattering length is much larger than the Fermi 
wave length but much smaller than the junction dimen- 
sions, the motion of the electrons in the normal layer is 
diffusive and the proximity effect may be described by 
the Usadel equations These equations are nonlinear, 
which complicates their analytical treatment except for 
several simple limits. One of the cases most accessible 
to an analytical treatment is the limit of a long disor- 
dered SNS junction (with the minigap energy scale much 
smaller than the superconducting gap) with transparent 
normal-metal-superconductor interfaces. The spectral 
properties of such a junction have been previously stud- 
ied in Refs. and we find it possible to further improve 
on those results. In this note we revisit this problem, re- 
fining some of the existing results and replacing numer- 
ical answers with analytical ones. We design this note 
as a quick reference on the structure of the minigap in 
a long disordered SNS junction which may be useful in 
view of renewed interest in such systems in connection 
with problems related to 7r-junctions ||^,^ and to meso- 
scopic fluctuations As a byproduct, we derive two 

useful identities for solutions to Usadel equations which 
simplify our analytical calculations. 

Assuming a quasi-one-dimensional geometry of the 
contact, the proximity effect in the normal layer may 
be described via the Usadel equations (in our paper we 
conform to the definitions of Ref . H ) 



^dle + ie sin6'- i(9^x)^sin26l = 0, 

a, (9,X sin^ 0) = , (I) 

where 0(x, e) and xi^^ ^-^e the variables parameterizing 
the zero angular momentum component of the Green's 
functions, g = cosO and / = sa\9 exp(ix); here, x is the 
phase of the superconducting correlations, and the local 
density of states p(x, e) (in units of the normal electron 
density in the bulk) is given by 

p{x,e) — Re cos9 . (2) 

Note that we measure lengths in units of the junction 
width L = 1 and the unit of energy e is the Thouless 
energy Ef. = Djl? , with D the diffusion constant in the 
normal metal. We further assume that the Thouless en- 
ergy i?c, as well as all other energy scales in the problem, 
are much smaller than the superconducting gap A. In 
this limit, the energy scale A is (to leading order) ex- 
cluded from the Usadel equations. 

Within the superconducting leads, Q = 7r/2 (at e ^ A) 
and X equals the phase of the superconducting order pa- 
rameter. Inside the junction, both Q and x turn into com- 
plex functions. For simplicity, we choose ideally transpar- 
ent interfaces and assume the normal metal to be much 
more disordered than the superconductor; in this case 
the boundary conditions become "rigid" @J^, 

e{x = W) = d{x = \) = ^-, 

X(x = 0) = 0; x(a; = l) = Xo. (3) 

The equations (|l|)-(||) form a closed set determining 
the density of states (with xo and e as input parame- 
ters). They have been analyzed in Refs. ||,^ and in this 
note we extend their results; we list them first and sketch 
their derivation afterwards. 

(i) The complex integral "density of states" defined as 

/5(e) = / dxcos0(x,e) (4) 
Jo 

obeys the relations 
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where B{e,xo) and C(e, xo) are integrals of (0), 
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(ii) Near Xo ^ and xo ~* ""j the phase dependence of 
the minigap Eg{xo) (in units of Ec) involves the leading 
terms, 



Egixo) = C2(l - Cixl) : Xo<7r, 
Egixo) ^ Csiir ~ Xo) , TT-xo^Ti-. 



(9) 
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The value of C2 has been derived in 



1J0 Jo (sinhiJo - sinh?9)i/2 



3.122, (11) 



with the maximum attained at i?o ~ 1.421. The value of 
Ci reported in |^ is incorrect by an order of magnitude 
and we find the correct value 
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0.0921 , 
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For C3 we find the analytic result close to the numerical 
value reported in 



C3 = — « 2.467 . 



(13) 



(iii) The integral density of states p{e) = Re p{e) has 
a square-root singularity at e = -Eg [|| , 



p(e) - C4(xo) 



1/2 



(14) 



with the coefficient C4 diverging at xo ^ tt as 

C4{xo) ^ a{n ~ xo)-'^' ■ (15) 

This asymptotic form has been found numerically in Refs. 

and we confirm it here analytically, together with 
the value for a. 
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2.494. 



(16) 



(iv) At energies s much higher than Ec but much lower 
than the superconducting gap, the leading corrections to 
the (integral) density of states are 



p(e, xo) « 1 - + e-^Ce (e) cos(xo) , (17) 



where 



C5 = 2 - 72 w 0.586 , 



Ce{e) = 16tan2 - 
8 



- cos ^/e ' 



cos ^/e — sin y/e 



(18) 



(19) 



The verification of (i) is straightforward: We integrate 
the Usadel equations (|l|) once and denote the coordinate- 
independent integrals by B and C as in (0) and (I). The 
relations (^ and (^) are then obtained via integration 
by parts in x and repeated use of the Usadel equations, 
see the Appendix for details. Relation ^ expresses the 
conservation of the total number of states as a function 
of Xo; indeed, (^ implies that p{e) is a total derivative in 
energy of the expression {2ep— iB) which is independent 
of Xo at large energies [this follows from our discussion of 
the result iv) below]; therefore, the integral of p{e) over 
energies is independent of xo- The identity (^ follows 
from the fact that the supercurrent (which equals ImC 
in appropriate units [^) can be expressed as the deriva- 
tive of the free energy with respect to xo; we verify (||) 
directly on the level of Usadel equation in the Appendix. 

To obtain the results of paragraph (ii) we follow the 
usual procedure in integrating the Usadel equations: be- 
low the gap, 6{x) takes the form — 7r/2 -I- id, where 
■d{x) is real and varies from at the N-S interface to 
its maximal value in the middle of the normal layer. 
From (0) and (||) we obtain the differential equations 



d^d^^fi-do)- f(d)/V~e, a,x = Ccosh-2^, (20) 

where /(i9) = sinh('i9) — (C^/4e) cosh~^ d. Their integra- 
tion over half the junction provides us with the solutions 
of the Usadel equations in the form 
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expressing s and xo in terms of the new parameters C /\/£ 
and do- The minigap Eg{xo) is defined as the maximal 
energy e c omp atible with this solution (with a real "^ix)). 
Equation ( |l2| ) is now easily obtained by expanding in C 
which is the small parameter near xo = 0, while the re- 
sult ( [1^ ) is obtained from an expansion in e, the small 
parameter close to xo = where the minigap vanishes 
(see also the derivation below). 

A more accurate expansion in small s is necessary to 
obtain the results of paragraph (iii). The density of 
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states near the gap edge (|TJ) is derived with the help of 
relation (^) . The key idea of the calculation is that most 
quantities are regular functions of the new parameters t^o 
and e/C^ [except for the point Xo = where the mini- 
gap vanishes]. The square- root singularity appears 
when inverting regular functions at the extremal point. 
Remarkably, in this way we can compute the density of 
states at the gap edge using the solutions to the Usadel 
equations below the gap and analytically continuing them 
to energies above the gap. 

Close to xo = TT the energy e is small and we can ex- 
pand in the small parameter S = e/C^, 

e = (tt^ cosh^ i?o) (5 - TT f i - ^\ e^'>"S^ + ... , (22) 
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We will show below that the solutions at the gap edge 
near Xo = ''^ involve a large parameter -do ^ 1, and hence 
we may keep only the leading terms in e^" in the coeffi- 
cients of S^. Next we invert (03) to find 



sinhiJo , 
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4 cosh^ -do 

and substitute the result into ( p^ ) to express the energy 
e as a power series in (tt — xa) a-nd as a function of i^q, 



-(tt - Xo)tanhi?o 
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(tt-xo)'. (25) 



The gap edge is given by maximizing e in ( pq ) as a func- 
tion of do; the corresponding value i?o maximizing e is 

Xc^(^~^) 84.08; (26) 
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we see that do is indeed large at the gap edge (for small 
TT — xo)j albeit only logarithmically and with a small pre- 
factor 1/3. Upon substitution into (Ea), this result also 
gives us the next-order correction to (|1C|), 



Egixa) = -r(7r- Xo) 
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We are now prepared to derive the density-of-states 



singularity at the gap edge ([14[)-(16). The maximum 
of the function £(190) at z9o is expressed in the relation 
e{i3o) ~ Eg — \d^g£\§g [t^o — '!9o]^/2; its inversion produces 
a square-root singularity at the gap edge in the function 



(28) 



and hence ?9o develops an imaginary part at energies 
£ > Eg. The singularity in i?o(e) translates into a square- 
root singularity in B{e); evaluating (^ in the junction 
middle and expressing C with the help of (p3) , we obtain 
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This singularity in B further translates, via (^, into a 
square-root singularity in the density of states p{e) . Eval- 
uating the coefficient in ( pof ) with the use of ( p5| ) for £(i9o) 
and of (H^) for iJoj we arrive at the final results (p^-(p^. 
The coefficient ( p^ agrees with the numerical findings in 

Another regime where the density of states is amenable 
to a simple analytic solution is at energies much larger 
than Ec but much smaller than A. In this limit, the cou- 
pling between the superconducting leads is weak, which 
allows us to derive the results of paragraph (iv). In equa- 
tion (p7|), the term proportional to is due to the 
suppression of the density of states near the interfaces, 
and the exponentially small term proportional to cos xo 
results from the Josephson coupling. 

At £ ^ 1, the Usadel equations may be solved by 
matching the solutions for the two semi-infinite N-S sys- 
tems; such a solution takes the form Q 



6'Ns(a;, £) = 4arctan tan , 



(30) 



where k = \J —2ie (Rck > 0) and the normal layer is 
at a; > 0. Matching of such solutions was performed in 
Ref. [|l^ by observing that for \Q\ <^ 1 (i.e., everywhere in 
the junction, except for the very thin layers at the inter- 
faces) the Usadel equations become linear in the variables 
/ = de^^ and / = 6*6^'^. Therefore the solution near the 
middle of the junction is given through the simple sum 

fix) = (x) + e-'^° 9^1 -x). (31) 
From this solution we easily find, using (Q) and (||), 

C(£) = 32tan2^Ke-'^sinxo, (32) 



B{e) 



-IE 



1 — 32 tan — e " cos xo 



(33) 



With the help of the identities (||) and (^) this imme- 
diately implies the results ( p^ ) and (|l^). The value of 
C5 is left undetermined by this method but may easily 
be obtained from directly integrating the local density of 
states (H) corresponding to the N-S solutions (30), which 
leads to the result ([l8|). 

We thank Urs Ledermann for discussions and for draw- 
ing our attention to Ref. |l^ and Swiss National Foun- 
dation for financial support. 
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Appendix: 

In order to verify Eq. (||) we use the definition of B to re-express 
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(34) 



The last integral contains five terms. After integrating the first term by parts in x (separating 89 /de), it vanishes 
against the third and fifth terms by virtue of the first Usadel equation (|l|). Integrating the fourth term in (^J) by 
parts in x annihilates it due to the second Usadel equation; we thus arrive at 



d_ 

8^ 



{iB-2ep) = -p, 



(35) 



which is equivalent to (H). 

In order to verify Eq. ^ we use the identity dx {8xX)'^ sin^ 9 — 2C + dx {8xX)^(^xo si^^ ^ to re-express 

8CI8e as 



.dC 



i 8 8 
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We use the second Usadel equation to differentiate the first of the two terms in e, 



(36) 
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In the right-hand side of Eq. (p7|), the first term is annihilated after integration by parts, and we finally arrive at 



^ir = -ol^ dxidxx)^sm29— + -— dx (dxx)^ sin29 — 
de 2 dxo Jo de 2 8e Ja dxo 
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where we have used the first Usadel equation. 



d9 
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= -2 



dx sin 9 



d9 ^ d ^ 
dxo 8x0 
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